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Abstract 

Multi species turbulence in inhomogeneous magnetised plasmas is found to exhibit symmetry 
breaking in the dynamical alignment of a third species with the fluctuating electron density and 
vorticity with respect to the magnetic field direction and the species' relative background gradi- 
ents. The possibility of truly chiral aggregation of charged molecules in magnetized space plasma 
turbulence is discussed. 
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I. INTRODUCTION: DRIFT WAVE TURBULENCE 

Drift wave turbulence Q is composed of coupled nonlinear fluctuations of pressure p 
and the electrostatic potential <fi in an inhomogeneous plasma with a guiding magnetic field 
B. The local electrostatic potential acts as a stream function for the dominant plasma 
flow with the drift velocity. The plasma is advected around the isocontours of a potential 
perturbation perpendicular to the magnetic field and forms a quasi two-dimensional vortex: 
an initial localized plasma density perturbation rapidly loses electrons along the magnetic 
field B, leaving a positive electrostatic potential perturbation <fr in its place. The associated 
electric field E = — V0 is directed outward from the perturbation centre, and in combination 
with a magnetic field B causes an "E-cross-B" vortex drift motion of the whole plasma along 
isocontours of the potential with a velocity = E x B/£? 2 . 

This gyration averaged drift motion is a result of the combined electric and magnetic 
forces acting on a charged particle with small gyration radius compared to background 
scales jj. In the presence of a background density (or pressure) gradient Vn the vortex 
propagates in the direction perpendicular to both Vn and B. Drift waves become unstable 
if the electron density can not rapidly adapt to the electrostatic potential along the direction 
parallel to B, which can be caused by interaction of the electrons with other particles or 
waves. Nonlinear coupling of drift waves results in a self-sustained fully developed turbulent 
state Q, ^. 

The requirements for the occurrence of drift waves are ubiquitously met in inhomogeneous 
magnetized plasmas when the drift scale p s = (MjT e ) / (eB) is much smaller than the 
background pressure gradient length L p = (Vlnp) -1 , or equivalently, when the fluctuation 
frequency u is much smaller than the ion gyration frequency u C i = eB/Mi. Here Mj is the 
ion mass, T e the electron temperature (in eV), e the elementary charge, B the magnetic field 
strength, and p = n e T e the electron pressure. 

Drift wave vortices are mainly excited in the size of a few drift scales p s and nonlinearly 
dually cascade to both smaller and larger scales. This perpendicular spatial scales of turbu- 
lent fluctuations are characteristically much smaller than any parallel fluctuation gradients 
and scales, so that wave vectors fulfill k\\ -C k±. 
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le aggregation and transport in turbulence is a widely studied subject in fluid dy- 
namics {4]. In turbulent plasmas the aggregation of charged particles is in addition influenced 
by static and dynamic electric and magnetic fields. For gyrofluid drift wave turbulence in 
inhomogeneous magnetized plasmas it has been observed by Scott [5] that the gyrocenter 
density of light positive trace ions (in addition to electrons and a main ion species) in three- 
species fusion plasma drift wave turbulence tend to dynamically align with the fluctuating 
electron density, so that fluctuation amplitudes of the electron and trace ion densities are 



spatiotemporally closely correlated. Priego et al. [10] have found that the particle density of 
inertial impurties in a fluid drift wave turbulence model is also closely correlated with the 
ExB flow vorticity. 

The present analysis shows that the dynamical alignment in drift wave turbulence is 
not only sign selective with respect to the vorticity of trapping eddies, for any given trace 
species charge and background magnetic field direction, but also their large-scale background 
gradient. 

The resulting asymmetric effects on transport and aggregation of charged particles is 
specifically discussed for application to space plasmas, in particular for cases of a nearly 
collisionless plasma with cold ions and low parallel electron resistivity, where the electrostatic 
potential closely follows changes in the electron density. It is shown that in this case drift 
wave turbulence constitutes a unique mechanism for truly chiral vortical aggregation of 
charged molecules in space environments with a background magnetic field and plasma 
density gradients. 

This paper is organized as follows: in section II the cold-ion gyrofluid model equations for 
the present multi-species drift wave turbulence studies are introduced. Section III describes 
the numerical methods, and in section IV the computational results are presented. Section V 
discusses implications of the results for chiral particle aggregation in turbulent space plasmas, 
which poses a possible extraterrestrial physical mechanism for achieving enantiomeric excess 
in prebiotic molecular synthesis. The relation between the present gyrofluid model and the 
inertial fluid model by Priego et al. [hJ is considered in the Appendix. 
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II. MODEL: MULTI-SPECIES (COLD ION) GYROFLUID EQUATIONS 



Multi-species drift wave turbulence in the weakly collisional limit can be conveniently 
treated within a gyre-fluid model. Here a strongly reduced variant of the gyrofluid electro- 
magnetic model ("GEM") by Scott [6( is adapted for numerical computation of three-species 
drift wave turbulence in a quasi-2D approximation in the electrostatic, isothermal, cold ion 
limit (neglecting finite Larmor radius effects) with a dissipative parallel coupling model. 

The gyrocenter particle densities n s for (electron, main ion, and trace ion) species s G 
(e,i,z) are evolved by nonlinear advection equations 

D t n s = (d t + V£ • V)n s = d t n s + [0, n s \ = C s , (1) 

where all C s include hyperviscous dissipation terms z/ 4 V 4 n s for numerical stability. The 
parallel gradient of the parallel fluid velocity components are for simplicity expressed by 
assuming a single parallel wave vector k\\ in the force balance equation along the magnetic 
field [2|. Then C e in addition includes the coupling term d(4> — n e ), where d is the parallel 
coupling coefficient proportional to k?Ji/ e i. For weakly dissipative plasmas d is well larger 
than unity, but for practical purposes d = 2 already sufficiently supports nearly adiabatic 
coupling between n e and while allowing reasonable time steps. The coupling can be 
regarded as non-adiabatic when d < 1. 

The electrostatic potential is here derived from the local (linear) polarization equation 

PmV 2 = n e - Hi - a z n z , (2) 

where p m = 1 + a z fi z with a z = Zn z /n e and fi z = m z / '(Zrrii). In the local model only the 
density perturbations are evolved in the advection equations, which then gain an additional 
background advection term g s d y 4> with g s = (L^/L ns ) on the right hand side, where L±_ is a 
normalising perpendicular scale (here set identical to L ne so that g e = = 1) and y is the 
coordinate perpendicular to both the magnetic field and the background gradient directions. 
In the following only trace ions with a 2 -C 1 are considered, and the parameter a z = 0.0001 
is going to be fixed, while other parameters are being varied. 

The 2D isothermal gyrofluid model and its relation to (inertial) fluid models is briefly 
discussed in the Appendix. 
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III. METHOD: NUMERICAL SIMULATION AND ANALYSIS 



The numerical scheme for solution of the advection and polarization equations is described 
(for a similar "Hasegawa Wakatani" type code where the ion density equation is replaced 
by a vorticity equation) in refs. jl, 8]. 

The present computations use a 512 x 512 grid corresponding to (64 p s ) 2 in units of the 
drift scale. This is a resolution well established for drift wave turbulence computations, 
resolving both the necessary small scales just below the gyro radius (or drift scale), and also 
allowing enough large scales for fully developed turbulent spectrum. The results converge 
for both higher resolution and same size in drift units, and for larger size with accordingly 
increased grid nodes. 

The computations are initialized with quasi-turbulent random density fields and run until 
a saturated state is achieved. The resulting turbulent advection dynamics is restrained to 
the 2-dimensional plane perpendicular to the (local) magnetic field [[J. Statistical analysis is 
performed on spatial and temporal fluctuation data in this saturated state. The turbulence 
characteristics of drift wave systems (like power spectra, probability distributions, etc.) in 
general have been extensively discussed before, so that here the focus is on analysis of the 
correlation of the fluctuating density of the additional species with the other fluctuating 
fields. 



IV. RESULTS: ASYMMETRIC DYNAMICAL ALIGNMENT 

In accordance with ref. jlj] we also observe in our simulations that the gyrocenter density 
of a low-density massive trace ion species tends to dynamically align with the electron density 
n e . 

The anomalous diffusion, clustering and pinch of impurities in plasma edge turbulence has 
also been studied previously by Priego et al. 10] . In this paper the impurities were modeled 
as a passive fluid advected by the electric and polarization drifts, while the ambient plasma 
turbulence was modeled using the 2D Hasegawa- Wakatani model. As a consequence of 
compressibility it has there also been found that the density of inertial impurities correlates 
with the vorticity of the ExB velocity. Trace impurities were observed to cluster in vortices 
of a precise orientation determined by the charge of the impurity particles 10 ], 
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The major difference to the present approach is that Priego et al. have considered only a 
constant impurity background (i.e. no impurity background gradient), but have in addition 
included additional nonlinear inertial effects through the impurity polarization drift. Linear 
inertial effects by polarization are already consistently included in the local gyrofluid model. 
The detailed correspondence between the two models is discussed in the Appendix. 

For an initial homogeneous distribution of impurities with no background gradient (g z = 
0) and vanishing (or constant) initial perturbation n z (t — 0) =0 the perturbed gyrocenter 
density n z (t) remains zero for all times, which directly follows from the evolution equation 
Dtn z = 0. The perturbed fluid particle impurity density N z is related to the cold impurity 
gyrocenter density by the transformation N z = n z + fi z Q (see Appendix). As a result the 



perturbed fluid particle impurity density evolves according to D t N z = [i z D t Q and direct 



y 



follows the vorticity, which has also been observed in fluid simulations by Priego et al. [10]. 
In the following the focus is on additional effects of a finite impurity background density 
gradient specified by g z . 

In the adiabatic (weakly collisional) limit the electrons can be assumed to follow a Boltz- 
mann distribution with 

(n + n e ) = n exp[e(p/T e ] « n (1 + e<p/T e ), (3) 

so that (for given T e and n ) a positive n e perturbation corresponds to a positive localized 
potential perturbation 0. The ~ (— x B) drift motion of a plasma vortex around a 
localized <fi perturbation possesses a vorticity 

= Vxv £ = (B/|S|)VV, (4) 

related to the Laplacian of the electrostatic potential, and thus (depending on the sign of the 
fluctuating potential) a definite sense of rotation with respect to the background magnetic 
field B/\B 



In ref. the absolute correlation \r(n e , n z )\ between the (gyrocenter) density pertur- 
bations of electrons n e and trace ions n z has been determined. Our present analysis in 
addition shows that under specific conditions a definite sign relation appears for the sample 
correlation coefficient 
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FIG. 1: Left: Vorticity Q(x,y); right: trace ion gyrocenter density perturbation n z (x,y). 
Normalized amplitudes are used and only a quarter of the complete computational domain 
is shown. Q and n z here show close negative alignment for weak (quasi- adiabatic) coupling 
with d = 2. The (co-aligned) impurity gradient length is here set to L±/L nz = +0.001. 



where the sum ^2 is taken over all grid points of the computational domain, and the bar n s 
denotes the domain average of the specific particle density. 

In local computations with the present model, where the densities are split into a static 
spatially slowly varying background component n Q with perpendicular gradient lengths L n = 
(Vlnno) -1 and a fluctuating part with small amplitudes n (typically in the range of a few 
percent), the sign and value of r(n e ,n z ) ~ ±(0.90 ± 0.02) are observed to only depend on 
the relative sign but not the magnitude of the gradient lengths L ne and L nz , for a given 
direction of the background magnetic field and fixed other parameters (d — 2, \i z — 10). 

The result for r(n e ,n z ) changes only marginally for most other parameter variations. 
In particular, the sign of the impurity charge Z has no effect on the alignment property. 
Stronger adiabaticity leaves r(d = 10) fa ±(0.90 ± 0.02) largely unchanged, while a smaller 
(strongly non-adiabatic) dissipative coupling coefficient results in r(d = 0.01) ~ ±(0.79 ± 
0.02). In Fig. [1] the computed turbulent fields of vorticity (left figure) and a heavy (for 
example molecular) charged particle species density perturbation (right) in the 2D plane 
perpendicular to B are shown in a snapshot to be closely (negatively) spatially correlated 
for parameters d = 2 (quasi-adiabatic) and g z = L±/L nz = ±0.001 (co-aligned background 
density gradients). 
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Table 1: correlation coefficient r(n e ,n z ) 
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The observed sign-selective multi-species dynamical alignment effect is basically caused 
by the (linear) drive of density fluctuations dtn z ~ g z d y <f) by this background advection 
term, where the potential <p is f° r each species just acting on the respective background 
gradients with length L ns . The species gyrocenter densities are enhanced (positive partial 
time derivative) when the background advection is positive, and decreased when negative. 
Co-aligned background gradients of electrons and trace ions thus dynamically also imply 
co-alignment of the gyrocenter density perturbations. In case of adiabatic electrons this 
additionally implies negative alignment with vorticity, as can be seen in Fig. [TJ For counter- 
aligned background density gradients the plot for n z (x,y) in Fig. [1] has exactly the same 
topology only with signs reversed, i.e. blue and red in the colour scale exchanged. 

Strong non-adiabaticity (d < 1) only slightly reduces the electron-to-impurity density 
correlation coefficient r(n e ,n z ) compared to adiabatic cases. Although the alignment be- 
tween electron density and electrostatic potential is reduced for a non-adiabatic response, 
the alignment between the impurity gyrocenter density n z and vorticity Q still remains. 

The results for r(n e , n z ) as a function of the dissipative coupling parameter d and (co- or 
counter-aligned) impurity gradient length g z = L±/L nz are summarized in Table 1. 



S 




FIG. 2: A Gaussian trace ion puff placed in the turbulent plasma is initially localized 
around the center of the x domain (and constant in y). Left: impurity density n z (x,y). 
Right: the product function fl(x, y) ■ n z (x, y) shows in comparison with the top figure that 
impurity density and vorticity are preferentially negatively aligned in the right half of the 
domain (where the initial density gradients are co-aligned) and positively aligned on the 

left. 




100 200 300 400 500 100 200 300 400 500 



FIG. 3: A Gaussian trace ion puff placed in the turbulent plasma is initially localized 
around the center of the y domain (and constant in x). Left: impurity density n z (x,y). 
Right: the product function Q(x,y) ■ n z (x, y) shows that in this case no preferential sign of 
alignment between impurity density and vorticity is found. 
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The situation is different if there is no background distribution of the trace ion species, 
but only a smaller localized cloud (diffusing over time) with an initial spatial extension 
in the same order of magnitude as the turbulence scales. Then the cloud can have global 
gradients of its density in all directions with respect to the background electron (and primary 
ion) gradient, and the signs of r approximately can cancel to zero by integration over the 
computational domain, while the absolute correlation coefficient |r| remains near unity. 

The turbulent spreading of such a initially Gaussian trace ion puff (with zero background) 
strongly localized in x or y is shown in Figs. [2] and [3j respectively. In Fig. [2] a Gaussian 
trace ion puff placed in the turbulent plasma is initially localized around the center of the x 
domain and constant in y. The product function Q(x,y) ■ n z (x,y) shows (in the left part of 
the figure) in comparison with the impurity gyrocenter density n z (x,y) (right figure) that 
n z (x,y) and Q(x,y) are preferentially negatively aligned in the right half of the domain, 
where the initial density gradients are co- aligned, and positively aligned on the left. This 
behaviour is in accordance with the previous result for a constant background gradient, as 
has been shown in Fig. [TJ For the x-symmetric case in Fig. |3]no preferential sign of alignment 
between impurity gyrocenter density and vorticity is found. 

Also this case the perturbed fluid particle density N z is still correlated with vorticity 
according to the relation D t N z = ft z D t Q. The preceding discussion of the present results 
was funded within the gyrocenter density n z representation of the gyrofluid model. It remains 
to be clarified when the advection by the background gradient or the inertial contribution is 
the dominant factor that determines alignment of the fluid particle density N z = n z + [x z Q 
with vorticity. 

In Fig. H]the correlation coefficient r(Q, N z ) is shown as a function of the impurity density 
gradient parameter g n at d = 2 for \i z = 10 (circles, bold solid line), \i z = 30 (squares, thin 
dashed line) and fi z = —10 (diamonds, thin solid line). The alignment of perturbed impurity 
particle density with vorticity changes sign around g n ~ /i z /2: for g n <C ji z the alignment 
directly follows the vorticity with r(Q,N z ) ~ 1, and for g n ^> fi z the alignment is reversed 
towards negative vorticity with r(Q,N z ) rj —1. 

The computations have been repeated for a non-adiabatic parallel coupling coefficient 
d = 0.1, for which the results remain very similar. The overall impression given by Fig. H] 
does not change signficantly in this case, except for slight (order of per cent) differences 
within fluctuation averaging errors. 
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FIG. 4: Correlation coefficient r(Q, N z ) as a function of the impurity density parameter g n 
for n z = 10 (circles, bold solid line), fi z = 30 (squares, thin dashed line) and u, z = —10 
(diamonds, thin solid line): the alignment of impurity particle density with vorticity 

changes sign around g n ~ [x z /2. 

The sign of alignment is thus determined by the strength of the impurity gradient in 
relation to the mass-to-charge ratio of the impurities. For example, for singly positively 
charged heavy molecules (like typical space biomolecules) with /i z >> 1 and small impurity 
gradient g z ~ 1 the alignment is tendentially towards positive vorticity. For vanishing 
impurity gradient (g z = 0) the correlation r(Q,N z ) is always exactly +1 for positively 
charged impurities, and —1 for negative impurities. 

The basic conclusion is that for given background gradients and magnetic field direction 
the perturbed trace ion species gyrocenter density dynamically aligns with a definite sign 
of the electron density fluctuations, and thus of the electrostatic potential fluctations, and 
consequently of vorticity Q: for example, an excess of trace ions aggregates within vortices 
of clockwise direction, and a deficit is found in vortices with anti-clockwise direction (or vice 
versa, depending on global parameters). 

While, as usual in a fully developed turbulent state, vortices of both signs appear equally 
likely and evenly distributed over all turbulent scales, the trace particle aggregation on drift 
scales emerges with one prefered rotationality with respect to the background magnetic field 
and background particle gradients. 
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V. CHIRAL MOLECULAR AGGREGATION IN DRIFT WAVE TURBULENCE 



Finally, a possible relevance of the asymmetric alignment effect on molecular chemistry 
in a magnetized space plasma environment is suggested. 

Homochirality of biomolecules - the fact that the essential chemical building blocks of life 
have a certain handedness while synthetic production leads to equal (racemic) distribution 
of left-handed and right-handed chiral structures - has ever since its discovery by Pasteur 



12j] posed a formidable puzzle [3, Q . Pasteur already had unsuccessfully tried to identify 



physical causes for this biological symmetry breaking, for example by imposing chirality 
through fluid vortices in a centrifuge, and by exposing chemical solutions to a magnetic field 



15|- 



A number of theories and experiments on the origin of chirality have since been put 



forward, like effects of circularly polarized light on the molecular reactions 16Ml8| or a 



combination of magnetic field and non-polarized light [19j, or possible electroweak effects 
on quantum chemistry 20j. All of these mechanisms could be active in interplanetary 
and interstellar space, and would hint on an extraterrestrial origin of early fundamental 
biomolecules. 

That chirality can in principle indeed be imposed by rotational forces has been confirmed 



experimentally for different situations 



21 



22j, but what mechanism could invoke a specific 



directionality in turbulent (terrestrial or space) fluids or plasmas, where vortices of both 
senses of rotation generically occur mixed across all scales, has been left an open question. 
Here we argue that rotation asymmetric ion aggregation in drift vortices in magnetized space 
plasmas constitutes a mechanism for fostering a truly chiral environment for enantiomeric 
selective extraterrestrial formation of biomolecules. 

The drift scaling conditions and possibilities for occurrence of drift wave turbulence are 
well fulfilled for a number of typical space plasma parameters and magnetic field strengths, 
like in (warm ionized) clouds in the interstellar medium 23]. 

The background pressure or density gradient length of molecular ion species in the inter- 
stellar medium can take values across many orders of magnitude. While interstellar clouds 
range in size between a few and hundreds of parsecs (1 pc = 3.0857 • 10 16 m) and thus have 
global gradient lengths of the same order, the local gradients can be set by macroscopic 
turbulence (resulting from large scale magnetohydrodynamic motion or by chaotic external 



12 



drive through winds and jets) and vary widely, from observable scales of 100 pc down to 



1000 km 



Then again, the gradient of molecular ion density can also be set by a spatially varying 
degree of ionization, for example through inhomogeneous irradiation at the edge of molecular 
clouds. Magnetic field strengths in the interstellar medium can occur up to a few pG to mG, 
and temperatures range between 10 K in cold molecular clouds to 10000 K in warm ionized 
interstellar medium. The resulting drift (and vortex) scale p s is in the order between a few 
hundred meters and a few hundred kilometers. 

Typical vortex life times are in the order between seconds and hours. The mean free 
path (parallel to a magnetic field) between particle collisions is in the order of 1000 km, 
so that the plasma is only weakly collisional. The drift ratio is in the maximum order of 
5 = p s /L p = 0.1 and smaller, thus well below unity. 

The relative fluctuation amplitudes (compared to background values of the plasma den- 
sity) in drift wave turbulence are typically in the order of the drift ratio and thus here in 
the range of a few percent or below. Compared to macroscopic flow-driven or magnetohy- 
drodynamic (MHD) turbulence, drift wave turbulence is most effective on smaller scales, in 
the order of p s , and with relatively small amplitude. Drift wave vortices in the size between 
a few hundred meters and hundreds of kilometers are thus expected to be present in most 
interstellar plasmas. 

Next we consider, if such vortices truly can account for a chiral physical effect. Chiral 
synthesis usually requires a truly chiral influence 25|, |26( : the parity transformation (P : 
x, v — >■ —x, —v) has to result in a mirror asymmetric state that needs to be different from 
simply a time reversal (T : t — > —t) plus a subsequent rotation by it. A purely 2D vortex can 
not exert a truly chiral influence, while 3D funnel-like fluid vortices in principle do. While 
a bulk rotation per se can not cause any direct polarizing effect on the reaction path for 
molecular synthesis [27], particle aggregation on a supramolecular level in vortex motion has 
been experimentally shown to be able to lead to chiral selection 21] . To date it however 
remains unclear how chirality from a supramolecular aggregate (clusters or dust) could be 
transfered to single molecules {28]. 

In the case of drift wave turbulence the parallel wave-like dynamics imposes a wave 
vector ku that under P changes direction with respect to the vorticity f2 and B (which 
themselves are pseudovectors and remain invariant under P), equivalently to a 3D vortex 
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tube. Drift wave turbulence is thus truly chiral, although it characteristically appears quasi- 
two-dimensional. In case of a large-scale density gradient Viin along B the (small) fc|| is given 
by this gradient length, otherwise it is determined by local fluctuations. Molecular ions are 
advected by the drift vortex motion, and are subjected to chiral aggregation onto neutral 
molecules, clusters or dust particles that are not participating in the plasma rotation. 

If in addition (as to be rather expected in most cases) a trace ion background density 
gradient in the direction parallel to B is present, the resulting parallel wave vector ku 
prescribes in combination with vorticity O a definite true chiral vortex effect. Chemical 
reactions in this system occur by collisions or aggregation of ions, that participate in the 
chiral vortex rotation, with neutral molecules or clusters of the quiescent (non rotating) 
background gas. The co- or counter alignment of B with O can further act as chiral catalyst. 

The chiral influence however changes sign in different sides or regions of the (large-scale) 
molecular cloud, when the relative direction of background particle gradient and magnetic 
field direction reverses. This implies that different regions of the interstellar medium favour 
different tendencies for chiral selective aggregation, and thus for a potential enantiomeric 
excess in molecular syntheses. The excess rate could be expected in the order of the relative 
density fluctuations in drift wave turbulence, of a few percent. 

A possible problem that could not be addressed with the present electrostatic model is 
whether chiral alignment may be broken by strong Alfvenic activity. Electromagnetic com- 



putations for moderate beta values (a few percent as for fusion edge plasmas) in ref. {5] for 
absolute correlation have however shown no significant deviation of the alignment character 
compared to electrostatic computations, so that chiral alignment may be expected to survive 
for finite beta. 

Chiral aggregation in drift vortices should therefore be feasible locally at least in lower 
beta regions of the interstellar or interplanetary media. The necessary conditions for chiral 
aggregation thus appear rather restrictive and cosmologically rare. It may be concluded 
that either the suggested mechanism for chiral synthesis is subdominant (compared to any 
other of proposed or yet unknown mechanisms), or that a chiral excess of molecules should 
cosmologically be a rather rare phenomenon itself. 
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Appendix 

The correspondence between the present cold-ion gyrofluid model with linear polarization 



10| including inertial nonlinear polarization effects is in 



and the HW model of Priego et al. 
the following briefly lined out. 

For a general discussion on nonlinear polarization and dissipative correspondence between 



29| by Scott. 



low-frequency fluid and gyrofluid equations we refer to ref. 

The Priego model (with notations and normalizations adopted to fit ours) consists of the 
2D Hasegawa-Wakatani (HW) model 

D t (N e -x) = d(<f)-N e ) (6) 
D t n = d(<j> - N e ) (7) 

where Q = V 2 is the vorticity, D t = d t + • V is the convective derivative with y# = 
(B/|S|) x V0, and D t x = —d y 4> in a local model where the background density gradient 
enters into the length scale normalisation. For simplicity we neglect in this presentation all 
dissipation terms (where Priego et al. have used normal second order dissipation, while we 
used fourth order hyperviscous dissipation terms). Here capital letters are used for the fluid 
particle densities N s for distinction to the gyrocenter densities n s of the gyrofluid model. 
The particle densities fulfill the direct quasi-neutrality condition N e = iVj for N z <C N i0 . 

The massive trace impurities are in Priego's model assumed to be passively advected, but 
due to their inertia respond to the velocity v 2 = + \ pz with the additional polarization 
drift velocity v pz = — /i 2 D 4 V0 in the global compressional impurity continuity equation 

d t N z + V • (iV 2 v 2 ) = (8) 

for the full impurity density N z (whereas only fluctuations of N e are evolved in the HW 
equation). For incompressional ExB velocity (when the magnetic field is homogeneous) this 
is equivalent to the full (global) impurity density equation used by Priego et al., which is 
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(in our notation) given in ref. 10] as 

D t N z - V • {(N z D t V<f)) = (9) 

with £ = fi z 5o. The factor 5q = p s /L ne stems from the normalization <\> — > 5 o _1 e0/T eO of the 
local HW model. 

This introduces a nonlinear polarization term into the dynamics of massive impurities. 
In a local approximation (corresponding to the assumption of small fluctuations on a large 
static background), where only the ExB advection is kept as the sole nonlinearity, this 
equation reduces to: 

D t N z = (N z0 D t n. (10) 

Now we derive a set of HW-like vorticity- density fluid equations from the (cold ion) 
gyrofluid model, and compare this with the Priego model. 

The present gyrofluid model consists of the local gyrocenter density equations 

D t (n e + g e x) = d(<p-n e ) (11) 
Dtim + gix) = (12) 
D t {n z + g z x) = (13) 

and the local quasi-neutral polarization equation 



■(r - \)<ft + a a T x n a 



0. 



(14) 



where r s = T s0 /T eQ , a s = Zn s /n e and /x s = m s / (Zrrii). The gyrofluid advective derivative 
D t n s = d t n s + w s ] includes the gyro-averaged potential ip — Ti0 with T = [1 + b]^ 1 



and Ti 



,1/2 



[1 + (1/2)6] 1 in Pade approximation, where b = — r s /i s V 2 . In Taylor 



approximation one would have To ~ [1 — b] and Ti ~ [1 — (1/2)6]. The species coefficients 
are r e — 1, a e — — 1, fi e — for electrons, and aj = 1 for ions. 

The background gradient terms g s = d x In n s0 fulfill (for zero background vorticity) the 
global quasi-neutrality condition g e = gi + a z gz ~ 9i when gi 1 for trace impurities. In 
the normalization of perpendicular length scales by L± = L n = |9 x lnn e o| _1 we then have 

9e = 9i = 1- 

The perturbed particle and gyrocenter densities are connected via the elements of the 
polarization equation as 29 j 

a s N s = dsTxris + fi s a s V 2 4> (15) 
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so that gyrofluid polarisation corresponds to the fluid particle quasi-neutrality condition 
£ s a s iV s = 0. 

From now on again cold ions with Tj = r z = are assumed, and thus if) = <p an d 
r = Ti — 1. The local polarization then reduces to p m V 2 = n e — rii — a z n z with p m = 
J2s a sf I s = 1 + a z p z . 

It can be immediately seen that the local cold-ion gyrofluid model already includes the 
(linear) polarization effect that has been added in the Priego model, when the gyrocenter 
density in the impurity continuity equation D t n z = is substituted by the fluid particle 
density: 

D t N z = D t a z n z + D t p z a z Vt = p z a z D t Vt (16) 

where a s = a s /Z = n s o/n e0 . This is identical to the linearized form eq.[10]of the Priego model 
(where the full impurity density is evolved so that on the right hand side fi z a z —> fi z S N z0 ). 
Completely passive trace ions would be achieved in the gyrofluid model for \i z = so that 

Pm = 1- 

A fully global set of nonlinear gyrofluid equations would include the nonlinear quasi- 
neutral polarization equation [3^], which (with M s = m s /mi) is given as: 

(M s n s V<f))+J2<lsFin s = 0. (17) 

s s 

This equation is expensive to solve numerically and requires e.g. the use of a multi-grid 
solver. Many codes therefore apply the linearized form, which can be treated by standard fast 
Poisson solvers. This is usually regarded as a valid approximation as long as the turbulent 
density fluctuations are small compared to the background plasma density. 

The trace of the (cold ion) global polarization eq. [17] again delivers the fluid-gyrofluid 
density relations: 

N s =n s + V- (M s n s V0). (18) 
From this we get the fluid impurity continuity equation 

D t N z = AV ■ (M 2 n,V0) = V • {M z n z D t V<f)) (19) 

where twice D t n z = has been used. The term on the right hand side however still contains 
the gyrofluid density n z . The approximation 

D t N z w V ■ (CN z D t V<f)) (20) 
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with the fluid particle density N z appears acceptable if either ( <C 1 or if only small turbu- 
lent fluctuations on a static background can be assumed, so that the additional gyrocenter 

'he resulting fluid 



io|. 



correction could be considered an order smaller than the other terms, 
impurity equation [2U] is then equivalent to eq. M as used by Priego et al. 

Summing up, we here retain the effects of a co- or counter-aligned impurity density back- 
ground gradient and linear polarization, but neglect nonlinear polarization effects. Priego 
et al. have considered a constant impurity background density, but due to their inclusion 
of a nonlinear polarization term were able to also treat additional nonlinear clustering and 
aggregation effects of impurities by inertia within vortices. 
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